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ABSTRACT 

Aims. We report on the redshift of the lensing galaxy and of the quasar QJ 0158-4325 and on the lens model of the system. 
Mettiods. A deep VLT/F0RS2 spectrum and HST/NICMOS-F160W images are deconvolved. From the images we derive the light 
profile of the lensing galaxy and an accurate relative astrometry for the system. In addition we measure the flux ratio between the 
quasar images in the Mg II emission line to constrain the mass model. 

Results. From the spectrum we measure the redshift of the lensing galaxy (z = 0.317 ± 0.001) and of the quasar (z = 1.294 ± 0.008). 
Using the flux ratio in the lens model allows to discard the SIE as a suitable approximation of the lens potential. On the contrary the 
truncated-PIEMD gives a good fit to the lens and leads to a time delay of At^_B=-14. 5+0.1 days, with Ho=73 km s"' Mpc"'. 
Conclusions. Using the flux ratio to constrain the mass model favors the truncated-PIEMD over the SIE, while ignoring this constraint 
leaves the choice open. 
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1. Introduction 

The study of gravitationally lensed quasars is relevant both to 
cosmology and to the determination of the total mass in (lens- 
ing) galaxies. Time delays between multiple quasar images de- 
pend both on the mass distribution of the lensing galaxy and on 
the Hubble parameter. Ho (Refsdal, 1964). The translation of an 
observed time delay into a value of Hq requires the modelling of 
the surface mass density of the lens. For that purpose, geomet- 
ric and photometric constraints are needed such as: (a) spectro- 
scopic redshifts for both the source and the lens components and 
(b) accurate relative positions of the lensed images and of the 
lensing galaxy. Additional constraints to be used are the flux ra- 
tios (FRs) between the quasar images. Indeed, the FRs are equal 
to the magnification ratios at the image positions, and hence de- 
pend on the mass distribution of the lens. However, the contri- 
bution of the lens (lensing galaxy and environment) must be dis- 
entangled from other sources of flux variation such as: differ- 
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ential galactic extinction, intrinsic quasar variability, microlens- 
ing and/or millilensing (de)magnification, all phenomena which 
have different time scales and wavelength dependencies (see e.g. 
Anguita et al. 2008 for more details). In the case of radio loud 
quasars this is rather straightforward, yet it does not apply to the 
whole quasar population. For optically selected quasar, the sepa- 
ration of the different sources of flux variation is more complex. 
In spite of this complexity, one must attempt it to derive values of 
the image FRs genuinely related to the strong lens. In the case of 
doubly imaged quasars (which constitute ~60% of the strongly 
lensed quasars known so far), it is particularly important to use 
the FR to constraint the lens mass model, as these systems are 
naturally the least constrained. 

In this paper we report on new observations of the lensed quasar 
QJ 0158-4325, and on a model of the lensing system. The quasar 
and the datasets relevant to our study are presented in §|2] In § 3 
the spectra and image deconvolutions are discussed. We mode 
the lens in § [4] Conclusions are given in § [5] 
We assume a WMAP type cosmology (Spergel et al. 2003, 
2007): //o=73 km s"' Mpc"', Q,„=0.3 and Qa=0.7. All mag- 
nitudes quoted in the paper are in the AB system. 

2. QJ 0158-4325 and dataset 

The quasar was detected in the Calan-Tololo Quasar survey 
(Maza et al. 1995), therefore called CTQ 414, and confirmed 
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to be a lensed quasar (at z=1.29) by Morgan et al. (1999). The 
two images of the quasar are separated by 1.22". A galaxy 
overdensity was detected along the line-of-sight to the quasar 
(Faure et al. 2004), suggesting the presence of a galaxy group. 
There remains in the literature a discrepancy of ~0.8" between 
the absolute declination of image A (see label in Fig. |3]l 
by Morgan et al (1999) and that provided in the CASTLES 
database (Munoz et al 1998). In order to elucidate this offset, 
we have matched the USNO-B 1 stellar catalogue to the FORS 
R-band images of the lens (presented in Faure et al. 2004): the 
RMS values are 0.07" and 0.06" in RA and DEC respectively. 
We have then derived the absolute astrometry of image A: 
ajim) =0l''58"'41.43' and (Jyzooo =-43° 25'3.5". 
Recently Morgan et al. (2008) have published the continuum 
light curves of the quasar images, monitored during four 
seasons at the Apache Point Observatory (SMARTS) and at 
ESO/La Silla (Euler Swiss Telescope). However, they were 
not able to derive a reliable time delay, as microlensing events 
overwhelmed the intrinsic variability, and their best guess, for 
a realistic mass model (image A leading) is: Af^^ — 20 to 
-10 days. The system is still under photometric monitoring. 
We have retrieved NICMOS-2 (Near Infrared Camera and 
Multi-Object Spectrometer) F160W images of QJ 0158-4325 
from the HST archival data related to a large public survey of 
gravitational lenses (PI: Falco, CASTLESQ). 
From the VLT archive, we have retrieved a set of four long-slit 
spectra centered on the quasar images and on the lensing 
galaxy, obtained with FORS2 (PI: Rix). Each spectrum has 
an exposure time of 1300 s. The spectra were taken through 
an 0.7"slit, under a mean seeing of 1.28". The resolution is 
4.86 Apix"', from 5500 A to 11000 A. The spectra have been 
reduced following Eigenbrod et al. (2006). It is not possible to 
disentangle the quasar light from the light of the lensing galaxy 
through a direct analysis of the spectra. However, the spectra 
are suitable for spectral deconvolution. 



3. Deconvolution 

3.1. Spectral analysis and flux ratio 

In order to measure the redshift of the lensing galaxy we have 
applied the spectral version of the MCS deconvolution algo- 
rithm (Magain, Courbin & Sohy 1998, Courbin 1999) to sep- 
arate the lensing galaxy from the much brighter quasar images 
in the FORS2 spectroscopic dataset. The long-slit FORS2 data 
do not sample any star To circumvent this difficulty, we have 
modified the deconvolution strategy, compared to previous stud- 
ies (e.g. Eigenbrod et al. 2008): we do not use stars, but instead 
the two bright quasar images as point spread function (PSF) ref- 
erences. We then deconvolve the data considering three point 
sources (two for the quasar images A and B, and one for the lens- 
ing galaxy). Along the deconvolution, the software optimizes 
the positions and intensities of the spectra, and at the end yields 
the individual spectra of the three point sources. The spectra of 
the quasar images display strong Mg II and weak [NeV] and H6 
emission lines (Fig.[T]). The bright Mg II doublet is used to derive 
the redshift, through a two-line fit. Using the measured line cen- 
ters and the rest wavelengths [2795.5A, 2802.7A] for the Mgll 
doublet, we deduce a redshift z=1.294±0.008 for the quasar 
Fig. |2] shows the final one-dimensional spectrum of the lens- 
ing galaxy, smoothed with a 24 A box. Residuals of the quasar 



Mg II line are negligible, assessing the efficiency of the spectral 
decomposition. The lensing galaxy spectrum exhibits the typical 
absorption lines of an elliptical galaxy at z=0.317±0.001 (based 
on the position of the following Unes: Gband, Hfi, Mgb, Fell, 
NaD and Ha). 
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Fig. 1. Spectra of the images of QJ 0158-4325. Top spectrum: 
image A. Bottom spectrum: image B. Bottom panel: ratio be- 
tween the two spectra. 



In §[T]we have listed the phenomena which can affect the ob- 
served FRs between quasar images. Let us discuss now the FR 
genuinely due to the lens and the confidence with which it can 
be used in the lens mass model. 

The quasar image spectra exhibit clearly the Mg II line, which 
is expected to arise from a region in the quasar, (i) which is too 
large to be sensibly magnified by microlensing (Abajas et al. 
2002, Eigenbrod et al. 2008), and (ii) with intrinsic flux varia- 
tion appreciable over a time scale of months to years (Peterson 
2001). Therefore, this broad, low ionization emission line is 
mainly affected by galactic extinction and by the lens magnifi- 
cation. Another source of flux variation could be substructures 
associated to the lensing galaxy on the milli-arcsecond scale 
(such as a dark matter clump or a dwarf galaxy). In the available 
dataset, we have no evidence for the presence of such substruc- 
tures. Therefore we assume that there is no perturbation of such 
scale. 

The FR in the Mg II line is a function of the flux measured in 
the emission line, and in the continuum, C,. The ratio be- 
tween images A and B spectra (Fig. [T] bottom panel) exhibits 
two remarkable facts, (i) a different value of the FR depend- 
ing on whether it is measured in the continuum or in the Mg II 
line and (ii) a slightly negative slope. As to the first fact, the 
quasar continuum is associated with the emission of an accretion 
disk, the projected size of which is comparable to the Einstein 
radii of stars in the lensing galaxy (e.g. Anguita et al 2008). 
Therefore, the continuum emission is prone to be magnified (or 
de-magnified) by individual stars in the lensing galaxy. The in- 
trinsic variation of this region may happens to be fast (days or 
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Table 1. Summary of astrometric, photometric and redshift 
^ ' ' ' ' ^ ' ' ' ^ measures. 



4 - H7 Hp Mgb Ha 




-I II I I l|l II I 

-|!! §' '4 B ' 

lAlj I I I I [L±l I I I U \ I I Ll I \ I I I L 

6000 7000 8000 9000 

X [A] 

Fig. 2. Spectrum of the lensing galaxy. The total integration 
time is 5000 s. The template spectrum of a redshifted elliptical 
galaxy at z=0.317 is shown for comparison (Kinney et al. 1996). 
Atmospheric absorptions are indicated by the label atm. 
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the de Vaucouleurs profile, as already mentioned in Rusin et al. 
(2003) and Morgan et al. (2008). We find the following parame- 
ters for the lensing galaxy light profile: PA =61°, e=0.32+0.01, 
re//=0.70±0.01". The magnitude and relative position are dis- 
played in Table [T] The value of the effective radius is in agree- 
ment with a previous study of this system by Rusin et al. 
(2003, Teff =0.66+0.04") for the lensing galaxy, using data from 
the HST WFPC (F814W and F555W) and the same NICMOS 
(F160W) observations. In the band F160W, they found the lens- 
ing galaxy to be ~0.2 mag brighter than in the present study: an 
explanation would be that the light of the host galaxy under the 
quasar image A was not removed properly, leading to an over- 
estimation of the lensing galaxy total flux. In Table [T] the error 
bars correspond to the standard deviations of the mean values 
measured in the individual deconvolved frames. 



weeks). As to the second fact, the slope in the continuum FR 
could be explained either by microlensing or by galactic extinc- 
tion. 

The continuum and emission line ratios vary for different rea- 
sons, hence to measure properly the FR in the emission line it 
is necessary to remove the contribution of the continuum. The 
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=5.52 (or AiniB 



FR in the Mg II line is given by: - - 
m^)=1.85 [mag]). The quasar image B being closest to the lens- 
ing galaxy center, it is at the same time less magnified by the 
lens and more extinguished by the lensing galaxy, than image A. 
Hence, the value A(otb - m^)=1.85 [mag] is an upper limit of 
the difference due to the potential of the lens. Yet, the elliptical 
nature of the lensing galaxy tells us that the differential galactic 
extinction is likely to be fairly small and that the lens magnifica- 
tion is indeed the dominant contributor. 



3.2. Image deconvolution and galaxy light profile 

Moreover, we have applied the MCS deconvolution code to the 
NICMOS/F160W images. A technique to build the PSF when 
no star is available in the field is the iterative method described 
in Chantry & Magain (2007). As a starting point, we used Tiny 
Tim PSFs (Kiist & Hook 1997) akeady improved by the itera- 
tive method applied on another lensed system (the Cloverleaf, 
see Chantry & Magain 2007) observed under similar conditions. 
Then we deconvolve simultaneously the four NICMOS expo- 
sures. Doing so we obtain a good estimate of the background 
level, which we subtract from the initial individual exposures. 
The images cleaned from their background are then deconvolved 
simultaneously to reach a final FWHM=0.07". The result is dis- 
played in Fig. [3] The quasar host galaxy is visible under quasar 
image A. 

In addition, a fit of the lensing galaxy luminosity profile is per- 
formed during the deconvolution. A good fit is obtained with 



4. Mass models 

We use the Lenstool code (Kneib et al. 1993, Julio et al. 2007) 
to model the lens, where the singular isothermal ellipsoid (SIE) 
is parameterized by its position, its position angle, its ellipticity 
and its velocity dispersion. During the modelling, the position, 
ellipticity and position angle of the lens are fixed to the values 



measured on the deconvolved NICMOS images (see § 3.2 1. The 



redshifts of the source and lensing galaxy are fixed to the new 
values measured in this paper. The only free parameter is the 
velocity dispersion. We find an acceptable fit (x^ ~18) for the 
SIE, but the predicted FR between the quasar images is four 
times larger than observed (see § |3.1| l. The predicted time delay 
is -12.2+0.1 days. If we force the FR to be equal to the value 
measured in the Mg II line, and in exchange let free the position 
angle of the galaxy or its ellipticity, we do not find any good fit 
(X^ >>100), showing that a single SIE cannot explain such a 
low FR between the quasar images. 

Then, we model the lens with a truncated pseudo-isothermal 
elliptical mass distribution (TPIEMD; Kassiola & Kovner 1993, 
Kneib et al. 1996) for both the bright and dark components of 
the lens. The TPIEMD has a core radius r^ore and cut radius r„„: 
we fix rcut^^Veff (see Eliasdottir et al. 2007) and rcore^O.Ol" , 
with the aim of keeping Xcore «^cut. By fixing r^,,, as a function 
of Xgff we effectively define a constant mass to light ratio lens 
potential. Doing so, the number of free parameters is reduced 
to one: the velocity dispersion. The FR predicted is slightly 
larger than the upper limit calculated from the spectra. If we 
force the FR to match the observed one and free the position 
angle or the ellipticity (within the error bars given by the 
deconvolution) we obtain a good fit {x^=2.\5). Using a different 
value for Xcore (equal to 0.005" or 0.05") modifies slightly the 
time delay (respectively -14.2±0. 1 days and -14.7±0.1 days). 
Results of the models are summarized in Table |2] The error bars 
correspond to Ay^^l. 
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Fig. 3. Deconvolution of QJ 0158-4325 images (F160W). North is to the top. East to the left, a) Original image, b) Deconvolved 
image: the quasar images A and B appear as point sources, the lensing galaxy is the sum of a de Vaucouleurs profile plus the 
residuals which are not fitted by the light profile. The white contours trace the light distribution of the smooth objects (quasar host 
galaxy and flux residuals of the lensing galaxy), (c) de Vaucouleurs luminosity profile of the lensing galaxy, with a lower flux cut 
than in image (b), in order to highlight its orientation and ellipticity. 



Table 2. Parameters of the lens mass models. Col. 2: of the 
fit. Col. 3: Velocity dispersion (km s~'). Col. 4: Time delay in 
days (image A leading). Col. 5: FR (image A brighter). 
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5. Conclusions 

We have measured the redshift of the lensing galaxy and im- 
proved the measurement of the quasar redshift in QJ 0158-4325. 
While modelling the lens, we noticed that using only the quasar 
image positions to constrain the mass model does not allow to 
discriminate a mass density profile over another: in the case of 
QJ 0158-4325, both a SIE and a TPIEMD appear as suitable ap- 
proximations of the lensing potential. If we add the constraint 
of the quasar image FR (measured in the Mg II line), it becomes 
possible to discard the SIE as a suitable approximation of the 
lensing potential. The lens modelled as a TPIEMD is good. The 
time delay inferred with the TPIEMD is AtA-B=-14.5±0.1 days 
for//o=73kms-i Mpc"'. 

In conclusion, using the FR to constrain the lens allows to select 
the most suitable galaxy profile and to measure its mass, as well 
as to predict the time delay with more confidence (or the Hubble 
constant if the time delay is known). 
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